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We study collective helping behavior and bystander effects in a coevolving helping network model. 
A node and a link of the network represents an agent who renders or receives help and a friendly 
relation between agents, respectively. A helping trial of an agent depends on relations with other 
involved agents and its result (success or failure) updates the relation between the helper and the 
recipient. We study the network link dynamics and its steady states analytically and numerically. 
The full phase diagram is presented with various kinds of active and inactive phases and the nature 
of phase transitions are explored. We find various interesting bystander effects, consistent with the 
field study results, of which the underlying mechanism is proposed. 
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I. INTRODUCTION 

In recent years, statistical mechanics has played an 
important role in quantitative investigation of the so- 
cial, economic, and psychological phenomena [Ij. Some 
social phenomena bear the resemblance to the physical 
processes in a sense that the macroscopic complex pat- 
terns emerge from the interaction of a large number of 
microscopic constituents. Along this line, one of the in- 
teresting subjects to physicists is opinion dynamics, most 
of which are concerned either with agent based model 
studies 0, Q or with the analysis of real data such as 
election results [1, [H]- For the model studies the effect 
of interaction structure has been considered in terms of 
the complex network theory Q. More recently, various 
coevolving networks with the agents interacting on those 
networks have been proposed and studied 0- The co- 
evolution of networks with agents may be justified for 
a broad range of social phenomena, because each agent 
may construct or reconstruct his or her own neighbor- 
hoods in response to the interaction with other agents. 

In this paper, we study a coevolving helping network 
model, known as the rescue model [1]. This model was 
motivated by the bystander effect observed in social psy- 
chology: The witnesses are less likely to intervene in the 
emergency situation when there are more witnesses. The 
bystander effect was introduced by Latane and his col- 
leagues in the late 1960's 0. Since then, the prosocial 
and helping behavior was investigated extensively from 
micro- to macro-level perspectives [lo| . The meso- level 
perspective corresponds to the study of helper-recipient 
dyads. We would like to see how the collective help- 
ing behavior at the level of society emerges from the re- 
peated helper-recipient interactions. We adopt the co- 
evolutionary dynamics consisting of two stages such as 
(a) an agent's behavior depending on the relations with 
other involved agents and (b) the relation update result- 
ing from the agent's behavior. 

The assumptions of the model are based on the results 
of social psychological experiments in the laboratory [§, 
[llj : First, the witnesses who had an acquaintance with 



the victim were faster to report the victim's distress than 
did the other witnesses. Second, the witnesses who had 
friends among themselves responded to the emergency 
situation faster. Thus the relations among the involved 
agents definitely affected their helping behaviors. For 
simplicity, we assume that the various relations among 
agents, such as friendship or short acquaintances, can be 
modeled by a simple (unweighted and undirected) link 
between nodes of a network, where a node represents an 
agent. 

We interpret the link as one node's expectation of help 
by the other node. The link density of a network can be 
interpreted as the aggregate expectation of help in a soci- 
ety. A new link between the intervening witness (helper) 
and the victim (recipient) may be created as a result 
of the successful intervention. Or the existing link be- 
tween them may be severed when the intervention fails, 
because the failure may reduce the future expectation 
of help between each other. Furthermore, to build the 
feasible model we adopt the cost-reward model for the 
witness's arousal to the emergency [l^- One of the basic 
assumptions is that the degree of arousal that a witness 
perceives is a monotonically increasing function of the 
perceived severity and clarity of the emergency. Here we 
will use the more explicit term, i.e. the degree of willing- 
ness to intervene, rather than the degree of arousal. 

In the previous works on the rescue model Q, the au- 
thors investigated the effect of the number of involved 
witnesses k on the collective helping behavior. They 
reported a non-monotonic variation of the link density 
(aggregate expectation of help) with k at some moderate 
values of the model parameters, which are partly con- 
sistent with the field study results 13|. In this work. 



we present the full phase diagram through the compre- 
hensive analytic study of the model with small k in the 
whole parameter space. We identify the various kinds 
of active and inactive phases and analyze the nature of 
phase transitions. Furthermore, we study the fiuctua- 
tions which are responsible for finite-size effects and cor- 
relations. Numerical simulation results for large k reveal 
various kinds (normal, reverse, complex) of bystander ef- 
fects depending on the parameter values. We propose the 
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underlying mechanism for these bystander effects, which 
is consistent with the numerical and analytical results. 

This paper is organized as following. In Sec. [TTl the 
coevolving helping network model is briefly introduced. 
In Sec. lIIIi the active and inactive phases and the transi- 
tions among them are identified by the analytic calcula- 
tions and confirmed by numerical simulations. In Sec. lIVl 
the various bystander effects are discussed. Finally, we 
summarize the results in Sec. |Vl 



II. MODEL 

We briefly introduce the coevolving helping network 
(CHN) model 18|. In the CHN, a node represents an 
agent and a link between two nodes represents a friendly 
relation between two agents. The links are unweighted 
and undirected, so the network is defined by the symmet- 
ric adjacency matrix {pij} where pij = 1 if two nodes i 
and j are connected by a link and otherwise. 

At each time step, an accident (emergency) occurs in- 
volving a randomly chosen agent v (victim) and also ran- 
domly chosen k agents (witnesses) from a population of 
N agents (nodes). The set of k agents for the victim v 
is denoted by A^. Each accident carries the degree of 
its severity represented by q, which is randomly drawn 
uniformly from [0,1). Each witness i d needs to over- 
come, at least, its nonnegative potential cost Ci to inter- 
vene in an accident (rescue attempt). Furthermore, the 
possibility of the rescue attempt by the witness should 
be enhanced by the friendly relation between the victim 
and the witness as well as the number of friends of the 
witness in the other k — 1 witnesses. 

Summing up all together, we may write the degree of 
willing of the witness i to intervene in the acci- 

dent occurred at time t on the victim v as 

Xiv{t) = q + ap^y{t) + [2py (i) - 1] - (I) 

where a > 0, b > 0, and A* denotes the set A^, excluding 
the witness i. The interaction term between witnesses en- 
hances if the witness i has more friends than strangers 
in A* , and vice versa. The magnitude of the interaction 
term increases with k as it can vary from b{k — f ) (all 
other witnesses are friends of the witness i) to —b{k — I) 
(all other witnesses are strangers to the witness i). So 
we expect that the interaction term may dominate the 
evolving dynamics over all other terms for sufficiently 
large k. 

We select the witness i ^ Ay with the largest value 
of the intervener [8] . If there are more than one 

witnesses with the largest XiD ^ OUG of them is selected 
randomly. The intervention attempt would succeed if 
Xiy > and fail otherwise, which is followed by the link 
(relation) update between the victim and the intervener. 
If successful, the friendly relation will be set up or main- 
tained between the victim and the intervener. Otherwise, 



they become or remain strangers to each other. The link 
update dynamics is summarized as 

p,y{t + M)^6[x,,{t)], (2) 

where 0{x) is a Heaviside step function with 6{x) = I 
for a; > and for x < 0. We choose the time step 
At = 1/L for convenience where L = ^^^""'"•^ is the total 
number of possible links between N nodes. At the next 
time step, the above procedure is repeated and the time 
t is incremented by 1 after L such iterations. 

We focus on the link density p{t) and the hole density 
u{t) as the order parameters, respectively: 

p(t) = I-ii(t) = i^p,,(i), (3) 

which should be proportional to the aggregate mean ex- 
pectation of help (intervention) in a society, which can 
be one of the important social features. 

Finally, in this work, we assume = c > for all i for 
simplicity, which may be valid for a sufficiently homoge- 
neous population. 

III. NONEQUILIBRIUM PHASE TRANSITIONS 
A. The case with a — and 6 = 

The a ~ b — case is special. Every witness's degree 
of willingness is identical q—c and the intervener 

is randomly selected among witnesses, regardless of the 
number of witnesses k. The probability for the successful 
intervention {xm > 0) should be equal to 1 — c for < c < 
I. Thus, probabilistically, the number of links increases 
by I with the rate W+ = (I — c) (I — p) and decreases by 1 
with the rate W- = cp. For c > 1, W+ = and W- = p. 

In the mean-field (MF) framework ignoring the 
stochastic temporal noise, the rate equation for the link 
density in the time-continuum limit is given as 

'i = W,-W., (4) 

where 

W+ = Fil^c)il-p), 

W. = F{c)p, (5) 

with F{x) = X for < a; < I, I for x > I, and for 
a; < 0. It is easy to show that p{t) behaves as 

p(i)=p(0)e-*-fp^(l-e-*), (6) 

with the steady-state density poo = hmt_j.oo p{t) as 

Poo = Fil ~ c), (7) 

which implies that there is a nonequilibrium phase tran- 
sition at c = I from an active phase (poo = f — c for 
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c < 1) into an inactive phase {poo = for c > 1). The or- 
der parameter exponent /3 is defined near the transition 
point in the active side as 



Po 



(8) 



with the reduced coupling constant e which measures the 
distance from the transition point. In this case, e — I — c 
and 13 ~ 1. 

In the inactive phase, the steady state is the vac- 
uum (no link state) where the dynamics becomes com- 
pletely dead, which is called as an absorbing (trapped) 
state. However, this transition is different from other 
well-known absorbing phase transitions [isj . since the 
vacuum is not absorbing in the active phase. Further- 
more, the system at the transition point (c — 1) is not 
critical, but absorbing: The link density decays exponen- 
tially {p{t) = p(0)e-*). 

One may derive the exact Langevin equation for large 
L, including the stochastic noise (in the Ito sense |16| ) 
through the usual Fokker-Planck formulation as 



dp 
'dt 



= fip) + V9ip)m, 



(9) 



where f{p) =W+-W- and g{p) = {W++W-)/L. ^(t) is 
a white noise with zero mean satisfying {£,{t)S,{t')) = 5{t— 
t') where (• • • ) is the noise average. The drift function 
f{p) = (1 — c) — p and the multiplicative factor g{p) = 
[(1 - c) + (2c - l)p]/L for c < 1, while f{p) = -p and 
g{p) = p/L for c > 1. In the L — oo limit, the noise term 
becomes negiligible and Eq. (|4]) becomes exact. 

The fluctuations due to the stochastic noise can be 
derived, using the Ito calculus [IB] such that 



d{G{p)) I .dG\ \ I d^G 



(10) 



where G is a general function of p. Then, the order pa- 
rameter fluctuation x defined as 



X = L[{p^)^ 
behaves in the steady state as 



Xo 



c(l 







c) for c < 1 
for c > 1 



(11) 



(12) 



Note that the fluctuation vanishes as the system ap- 
proaches the transition point (c = 1). The fluctuation 
exponent 7' defined as x ^ e^'^ is given as 7' = — 1. 



B. The case with a > and 6 = 

For a > and 6 = 0, we have the same rate equation 
as in Eq. @ with the link creation and annihilation rates 
as 



w+ = F(l-c)(l-p)^ 

= F(c-a)[l-(l-p)'=], 



(13) 




a 



FIG. 1: The phase diagram for the case with b — 0. The pa- 
rameter space (a, c) is divided into four regions by two tran- 
sition lines of c = 1 and c — a. 
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FIG. 2: (Color online) The steady-state link density poo along 
the a = 0.5 line for k — 2,4, and 6. 



respectively. Note that the a = O''" limit is singular except 
the case of fc = 1. 

It is straightforward to analyze the rate equation, 
which yields the steady-state density and its long-time 
dynamics. The parameter space (a,c) is divided into 
four different regions separated by two transition lines 
of c = 1 and c = a (Fig. [1]). In region I (a < c < 1), we 
have a finite link density as 



l-[(c-a)/(l-a)]i/'^ 



(14) 



(active phase). In region II (c < a and c < 1), we find 
the fully connected network; poo = 1 {paradise phase). 
In region III (c > a and c > 1), the network becomes 
completely disconnected; poo — (isolation phase). In 
region IV (1 < c < a), both transition rates W+ = W- = 
and no dynamics occurs; p(t) — p(0). Of course, the 
region IV is unphysical, so we focus on the three other 
regions and the transition lines between them. 

Approaching the I-III transition line (c = 1~) inside 
the active phase, the order parameter scales as 



(15) 



with Ai = 1/[A:(1 - a)], ei = 1 - c, and /3i = 1. Near the 
I-II transition line (c = a+), the hole density u = 1 — p 
plays a role of the order parameter, which scales as 



A 



J2 

2^2 



(16) 
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with A2 = (1 — a)~^/'^, £2 = c — a, and (32 = l/k. 

Note that the steady-state hnk density pao becomes 
smaller as k increases. In fact, one can easily show from 
Eq. that this holds for the whole region of the active 
phase, see Fig.[2l It implies that the society becomes less 
friendly (lower link density) with large k (more crowded 
society like big cities) in the long-time limit {normal by- 
stander effect). Therefore, one may regard the normal 
bystander effect as the result of iterative intervention at- 
tempts at the level of society. 

From the dynamic rules (transition rates) of our model, 
Eq. p3l) . it is easy to understand how this happens. A 
new link may be created only when there are no existing 
links between a victim and any of k witnesses, which we 
call the k-hole constraint. (Note that this applies only 
when a > 0.) The probability of finding this situation is 
proportional to (1 — p)'^, which monotonically decreases 
with k. So one can expect that the steady-state link 
density should be smaller with larger k. However, with 
interactions between witnesses (6 7^ 0), the fc-hole con- 
straint is not absolutely necessary to create a new link, 
which will be discussed later. 

Furthermore, suppose we are allowed to control c, for 
example, decrease c by giving all agents the same in- 
centive for intervention. The efficiency of our policy 
might be measured by the slope of the link density, i.e. 
e = — . As can be easily seen in Fig. [21 the policy be- 
comes quite effective near the paradise phase and more 
effective for larger k. However, it will be less efficient 
near the isolation phase and much lesser for larger k. 
These results imply that the incentive policy would not 
work well in unfriendly and crowded communities like big 
and rapidly developed cities, but may work efficiently in 
friendly and crowded communities like a guild. 

Now, we study the long-time dynamics in various re- 
gions. In the regions I and HI, the link density decays 
exponentially: p{t) — p^o ^ e^*/"^ with the characteris- 
tic time r such that t^^ = fc(l — a)^/'^(c — a)^~^/'^ in 
the active phase and r^^ = k{c — a) for c < a + 1 or fc 
for c > a -|- 1, respectively in the isolation phase. Along 
the I-III transition line, the link density still decays ex- 
ponentially with — k{l ^ a), which implies that this 
transition line is not critical but absorbing, i.e. the c — 1 
line belongs to the isolation phase. Approaching the I-II 
transition line (c — a+), the relaxation time diverges as 
T ^ £2 with the relaxation time exponent i^y = 1 — 1/fc 
for fc > 1. So the incentive policy works efficiently in this 
region, but its effect will show up quite slowly. 

In the region II (paradise phase), the long-time decay 
dynamics follows a much slower power-law; the hole den- 
sity u[t) ^ t^^ with the decay exponent (5 = l/(fc — 1) 
for fc > 1 except that u{t) ~ e"^^"'''* for fc = 1. With 
larger fc, the dynamics becomes extremely slower. The 
I-II transition line belongs to the paradise phase, mean- 
ing that the system exhibits the same power-law decay 
dynamics. 

As in Eq. ^ , one may write the exact Langevin equa- 
tion including the stochastic noise with the transition 



rates given by Eq. ([T3)) . Near the I-III transition line, 
one can expand in powers of p and obtain that J{p) = 
(l-c)-(l-a)fcpand g{p) = [{I- c) + {2c- a- l)kp]/ L 
in the active phase, while f{p) — — (c — a)kp and g{p) — 
(c — a)kp/L (c < a -I- 1) or kp/L {c > a + 1) in the 
isolation phase. Then, the steady-state fluctuations are 
analytically obtained, similar to the last subsection, as 

Xoo ^ Biep^ with Bi = (c- a)/[fc(l - a)^] and 7^ = -1 
in the active side and Xoo = in the isolation side. 

Near the I-II transition line, we need to write the 
Langevin equation in terms of the hole density u as 

^ = /» + VIRe(t), (17) 

where f{u) = {c—a) — {l~-a)u^ and g{u) = [(c — a) + (l + 
a — 2c)u^]/L in the active phase, and f{u) = —(1 — c)u'' 
and g{u) = (1 — c)u^ / L in the paradise phase. 

One can derive the fc-th order fluctuation analytically, 
using the Ito calculus of Eq. (fTO|) as 

^ L - {u)] = S^er^'" , (18) 

where B2 = (1 — c)/(l — a)^ and 72^*^'' = —1. From 
the simple scaling theory, one can relate the ordinary 
mean-square fluctuation Xoo — Xoo"* with x^<x^ through 
the exponent relation of 72 = 72*''''' -I- (fc — l)/?2. Hence, 
we find 72 — —l/k with /32 = 1/fc. 

Finally, we discuss the finite-size effects near the I- 
II transition line. The standard finite-size-scaling (FSS) 
theory predicts 

X« = £2"'"V('=' (e^LV-^^) , (19) 

where the scaling function behaves as (p^'^^x) ~ 0(1) for 

large x and <t){x)^^'^ ~ x''^ ' for small x. Fig. [3][a) shows 
an excellent collapse of various numerical finite-size data 
for Xxl at fc = 2 and a = 1/2 onto a scaling curve t/)*-^^ {x), 
with the choice of the FSS exponent j>2 = 1. The same 
is true for the ordinary fluctuation Xoo with the scaling 
function 4>{x) which behaves as 4>{x) ~ 0{1) for large x 
and 4>{x) ^ x'^^ for small x (Fig[3jb)). It follows that the 
order parameter u should scale as u ^ 1^-^2/1^2 — i^i/^ 
at £2 = 0"*". Note that the FSS exponent i?2 does not 
satisfy the standard hyperscaling relation as 72 -I- 2/32 = 

In the paradise phase including the I-II boundary, only 
the link creation process is allowed (no link annihila- 
tion dynamics). So there can be inflnitely many absorb- 
ing states where each node has some unconnected links 
(holes) and the number of those holes is less than fc for 
all nodes. Then, as there are no configurations satisfying 
the fc-hole constraint at all nodes, a link creation attempt 
is impossible and the system is frozen dynamically. The 
hole density in such absorbing states should be propor- 
tional to kN/L or equivalently u ^ L~^^^ in this region 
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FIG. 3: (Color online) The scaling collapse of numerical data 
for the fluctuations (a) a-nd (b) Xoo at A: = 2 and a = 0.5 
in the active phase (e2 > 0). 



(^2 ^ 0), "which is also confirmed numerically. In addi- 
tion, there is a discontinuous jump in the hole density for 
finite systems across the I-II transition line, due to the 
presence of this frozen dynamics in the paradise phase. 



C. The case with & > 

We no"w consider more realistic cases "with the interac- 
tion term {b > 0). In this case, each -witness's degree of 
"willingness in Eq. ([Ij is enhanced if he has more friends 
than strangers among the other witnesses and vice versa. 
Then, one may naively expect that the interaction drives 
a high link-density net"work to become a higher one, and 
a lo"w link-density net"work to become a lo"wer one. As a 
result, both the isolation phase and the paradise phase 
"would expand into the active phase. This is true at 
a — and, for sufficiently large b (> l/[2(fc — 1)]), the ac- 
tive phase squeezes do"wn to disappear (see, for example. 
Fig. 15). 

Ho"wever, for a > 0, the fc-hole constraint is in effect, 
"which suppresses the link density increase in general. Its 
effect is particularly big in the high density net"works, so 
the paradise phase may shrink. At the same time, the 
interaction term may loosen the constraint for large k 
and thus the naive picture can be restored. In fact, the 
fc-hole constraint becomes loose for a < 2b{k — 2), "when 
one of the "witnesses "without a link to the victim may 




1 b 

FIG. 4: The phase diagram for the case with a — and k — 2. 



take a higher degree of "willingness than one "with a link 
to the victim. So, a link creation dynamics is possible for 
a configuration not satisfying the fc-hole constraint. 

At fc = 2, the fc-hole constraint remains intact for any 
value of b and a > 0. The link density decreases even in 
the high link-density net"works. Then, the active phase 
expands into the paradise phase and also extends to the 
large a region. In contrast, for fc > 3 and small a (< 
26(fc — 2)), we find that the naive expectation holds and 
the active phase shrinks due to the expansion of both the 
isolation and the paradise phase. 

As discussed in the last subsection for 6 = 0, the fc- 
hole constraint is the key factor for the normal bystander 
effect: the steady-state link density poo decreases mono- 
tonically "with fc. For b > 0, ho"wever, as fc increases, the 
constraint becomes looser and the interaction term be- 
comes stronger (~ 5(fc — l)(2p — 1)) in the high density 
regime (p > 1/2). In this case, we expect the density 
increase "with fc for sufficiently large fc after the ordinary 
density decrease "with fc for small fc. These non-monotonic 
bystander effects "will be discussed more in the next sec- 
tion. 

No"w, "we consider the k ~ 2 and fc > 3 case in more 
details. 



k 



First, consider the a = case. One may easily sho"w 
that the transition rates are given as 

W+ = {l-p)[F{l-c-b){l-p) + Fil-c + b)p], 
= p[F{c + b){l- p) + F{c^b)p], (20) 

respectively. We find the phase diagram as in Fig. |4] 
"where the four different phases are separated by t"wo tran- 
sition lines of c = 6 and c = 1 — 6. As expected, the active 
phase shrinks "with increasing b by invasion of both the 
isolation and the paradise phase. In the active phase, 
the steady-state hnk density is Poo = (1 ^ ^ ^ c)/(l — 26), 
"which scales linearly (/3i = 1) near the Till transition line 
(c = 1 — 6) and the hole density u^o also scales linearly 
(/32 = 1) near the TII transition line (c = 6). 
With finite a > 0, the fc-hole constraint is in effect and 
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FIG. 5; The phase diagram for the case with k = 2 and < 
6 < I . The parameter space (a, c) is divided by four transition 
lines such as c — 1 ikb, c = a — b, and c = (1 + 2a — 6) /3. The 
dashed hues are 7 non-analytic lines for the transition rates. 



the transition rates are modified as 



W+ 



(1 - ~c~b)il-p) + F{l-c + b)p] , 

p{2-p) [F{c-a + b){l-p) + F{c-a-b)p], 

(21) 



respectively. Note that the a = 0+ hmit is again singular 
due to the fc-hole constraint. 

The parameter space (a, c) is divided into 14 different 
regions by 7 lines such as c = 1±6, c = b, c — a+(l±6), 
and c — a ± 5, where the transition rates change abruptly 
due to the non-analyticity of the function F(x) at a; = 
and 1. The divided regions are grouped into the four 
phases as shown in Fig. [SJ except that one central region 
is divided into two phases of I and III for < 6 < 1/2. 
As expected, the link density decreases with b for most 
regions [l7j , mainly due to the fc-hole constraint. The 
isolation phase III invades the active phase I and the 
paradise phase II retreats. We find the similar result for 
6 > 1/2 (not shown here). 

The active phase I comprises of four sub-phases where 
the order parameter p^ behaves differently: 

1 - Poo = 



' fc-|-^6^ + (l-a)(c-a-&) 

c—a-\-b 
l+a-2c 



l-Vl+4a'^ -4ab-4ac 
2a 



for c> a + b (la) 
for c > 1 - 6 (lb) 
for c < b (Ic) 

l+h-c-.,/4a2-362+2b+l-2(4a+b+l)c+5c2 . ,^ 

^ Tn — rr — ^ — ^ otlicrwise (Id 

2(a+o— c) ^ ' 

(22) 



All transitions are continuous but not difFerentiable be- 
tween the phases of I, II, and III, as well as between the 
four active sub-phases. In the vicinity of all transition 
lines, either the link density poo or the hole density Uoo 
vanishes linearly (/3i = /32 = 1). 

It is interesting to note the emergence of a nontriv- 
ial phase transition between lb and III in the region of 
1 — 6<c< 1-1-6, where the link creation and the link 
annihilation dynamics compete each other. In this re- 
gion, only the process representing the second term of 



W+ in Eq. ([2T|) is possible in the link creation dynam- 
ics. Therefore, the link creation dynamics is impossible 
when p becomes zero (no- link state). Once the system 
gets into the no-link state , it cannot escape out of that 
state {absorbing state). Based on these observations, we 
can niap the dynamics onto the well-known contact pro- 
cess [l8|, which exhibits an absorbing phase transition 
from vacuum (isolation phase) into an active phase. This 
transition is known to belong to the so-called directed 
percolation (DP) universality class [l5| . 

Near the transition between lb and III, we can easily 
write the Langevin equation for small p as 



dp „ „ 2 

- = ie,p-Bp 



ABp 



(23) 



where ei — c^^c with the transition point c^ — (l-|-2a — 
6)/3 and B = (1 — a + 26)/3 > 0. This equation is identi- 
cal to the MF Langevin equation describing the DP-type 
absorbing phase transition, which is characterized by the 
noise amplitude proportional to y/p. The steady-state 

link density behaves as poo — (3/i?)£i^ with (3i = 1. It is 
well known that the fluctuation exponent = in the 



active side for the MF DP universality class [15|. Uti- 
lizing the hyperscaling relation, we flnd the FSS expo- 
nent Di = 2, and thus expect poo ~ L^^i/'^i — L^^f^ 
at the transition [l^, which is confirmed numerically 
(not shown here). The relaxation time also diverges as 

T ^ e-^ with i^ii =1, in contrast to the case of the la-III 
transition line where r is finite. 

The transition from lb to II is also interesting. The cor- 
responding Langevin equation in the active side is given, 
in terms of the hole density u = 1 — p, as 



du 



iBu^ 



6Bu' 



(24) 



where £2=0 
hole density scales as u. 



2 with C2 = a — b. 



The steady-state 

P2 



(l/(3B))e^= with /32 



1. 

This type of the multiplicative noise linear in the activity 
field u has been studied extensively in literatures to de- 
scribe various physical systems including nonequilibrium 
wetting and synchronization phenomena for spatially ex- 
tended systems ,2Qj. Naive power counting which is ex- 
pected to hold in the MF systems yields the FSS expo- 
•nent 1/2 = 1 and 72 = —1. Thus, in the limit of £2 = 0"*", 
we have u ^ L^^. As u ^ L^^/"^ in the paradise phase 
and at the transition (£2 < 0), we expect a discontinuous 
jump in the hole density for finite systems at the Ib-II 
transition. Similar behaviors are found near the transi- 
tions from Id and Ic to II. All these results are confirmed 
numerically (not shown here). 

Note that the exponent /?2 changes from 1/fc = 1/2 
to 1 with the interaction term [b > 0). The relaxation 
time exponent i'^^ also changes from 1 — 1/A: = 1/2 to 
1. Inside the paradise phase II, the linear term (the first 
term in the right-hand side of Eq. p4))') vanishes and the 
hole density decay exponent becomes 5—1, which is the 
same as S = l/(fc — 1) = 1 at 6 = 0. The FSS exponent 
j>2 = 1 for both 6 = and 6 > 0. 
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FIG. 6: The phase diagram for the case with k — 3 and 



2. k>3 

For large k, we expect more complicated phase dia- 
grams. For example, see Fig. [6] for the phase diagram at 
fc = 3 and small b. The parameter space is divided by 11 
lines, such as c = 1 ± 26 , c = 1, c = 26, c = a + (1 ± 26), 
c = a + 1, c = a ± 2b, c = a, and a = 26. The last line 
of a = 26 is the boundary line between regions where the 
k-hole constraint is strictly valid and becomes loose. The 
divided regions are grouped into five phases including one 
new phase (inactive 'mixed' phase) in the central region 
(Fig.©. 

In this mixed phase, the system reaches either the iso- 
lation phase or the paradise phase, depending on the 
initial conditions and the stochastic dynamics. This is 
similar to the absorbing phase belonging to the directed 
Ising universality class 21] with two symmetric absorbing 
states. But, here, there is no symmetry between absorb- 
ing states. 

Our main interest lies in the active phase and the sur- 
rounding inactive (isolation and paradise) phases. As ex- 
pected, we find that the active phase shrinks due to the 
expansion of both inactive phases. But, the expansion 
of the paradise phase is both qualitatively and quantita- 
tively different in the regions of a > 26 and a < 26 at 
k = S. For a > 26, the k-hole constraint remains intact. 
The density decrease due to the k-hole constraint hap- 
pens to be balanced exactly by the density increase due 
to the interactions at A: = 3, which results in maintaining 
the I-II phase boundary (c = a) as it is at 6 = 0. For 
A: > 3, the interaction term dominates and thus the par- 
adise phase invades into the active phase. For a < 26 
where the fc-hole constraint is loose, the active phase 
shrinks more and the I-II phase boundary is given by 
the c = a -f 26 line. 

The hole density in the phase II decays as u{t) ^ t^^ 
with i5 = 1/2 for a > 26 and (5 = 1 for a < 26. For general 
k, we find i5 = l/[fc/2] {[x] is the integer value of x) for 
a < 26. Extending to the other regions of 26(n — 1) < 
a < 2bn with an integer n < fc — 2 is straightforward 
(not reported here). For a > 2b{k — 2) where the fc-hole 
constraint is valid, we find S = l/(fc — 1) (the same as 
that at 6 = 0). 



Finally, we report the values of scaling exponents near 
the I-II boundary. We find (32 — I, i>2 — l^, J^w = 2 for 
a < 26 and /?2 = 1/2, D2 = 1/2, = 1/2 for a > 26. 
Generalization to higher k is straightforward. Near the 
I-III boundary, we find /?i = 1, which holds for any k. 



IV. BYSTANDER EFFECTS AND SOCIAL 
IMPLICATIONS 

We investigate various bystander effects at the level of 
society/community for general fc, in relation to the em- 
pirical data analysis by social psychologists. Our model 
is set up, based on the experimental results in the labora- 
tory, but the specific values of control parameters cannot 
be directly inferred from these results. We compare our 
model study with corresponding field studies [ij, [l^l , in 
order to understand its social implications properly. 

The control parameters (a, 6, c, fc) can be regarded as 
given or intrinsic to a society /community, hence depend- 
ing on various factors of its socioeconomic condition and 
also on the taxonomy of helping behavior as well [13, Il4| . 
For the taxonomy of helping, social psychologists in- 
troduced three independent dimensions: doing (direct 
help) versus giving (indirect help), spontaneous versus 
planned, and serious versus nonserious. Helping mea- 
sures, such as helping rate, were observed from field stud- 
ies performed in the streets of 55 cities and towns in Aus- 
traha ^ and of 36 cities in U.S. [14]. 

For example, in one of experiments the experimenter 
walking on a street drops a pen accidentally in front of 
a solitary pedestrian, and continues walking, then sees 
whether the pedestrian picks the pen and bring it to the 
experimenter or simply ignores. The helping rate may 
correspond to the success rate in our model, which is de- 
fined as the number of successful interventions divided by 
the total number of interventions. In general the success 
rate is a monotonic function of the link density p, i.e. the 
aggregate expectation of help in a community. Therefore 
we can interpret p as a helping measure of a commu- 
nity. The helping behavior in this experiment belongs 
to the direct, spontaneous, and nonserious type with a 
low potential cost (small c) , while the rescue attempts in 
emergency situations are classified into the direct, spon- 
taneous, and serious type with large c. 

The socioeconomic factors {community variables) in- 
clude the population size, population density, cost of liv- 
ing, unemployment rates, and etc. Here the number of 
witnesses fc per accident corresponds to the population 
density (crowdedness) and the system size N to the pop- 
ulation size, respectively. While the population density 
and the population size are highly correlated in reality 
(big cities are usually crowded), we considered them as 
being independent to each other for generality. 

The previous empirical studies showed that the most 
relevant factors for the helping rate are the population 
density and the population size. The strong negative cor- 
relation between the helping rate and the population den- 
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FIG. 7: (Color online) The various bystander effects: (a) NBE, (b) CBE, and (c) RBE when b = 0.05 and a = 0, b, and 36. 
(d) Critical potential cost lines CR(a) and CN(a) in the (a, c) phase diagram at fixed b = 0.05. We find the NBE for c > CN(a) 
and the RBE for c < cr. In between these two boundaries, the CBE is found. When identifying the BE's we consider only the 
active phase, where only one stationary link density exists. 



sity/size was found. This implies the normal bystander 
effect, i.e. p decreasing with k in our model. However, 
in some nonserious experiments, the helping rate first 
decreases and then increases as the population size (or 
density) increases fisj, and a slightly positive correla- 
tion between the helping rate and the population density 
was also found 14]. These unusual bystander effects can 
be supported by various behaviors of p observed in our 
model study (see below). 

The experiments in the field studies have been done for 
measuring the aggregate expectation of help in a given 
society, not for simulating our model dynamics. We as- 
sume that the helping network in the society has evolved 
by our model dynamics and has become already station- 
ary at the time of the experiments. By comparing the ex- 
perimental data with the stationary results of our model, 
one may check the validity of our coevolving helping net- 
work (CHN) model and also predict some social features 
in real social systems. Especially, we focus on the various 
bystander effects found in the field studies. 

Now we discuss the fc-dependent behavior of the link 
density p, based on the numerical results for the CHN 
model, shown in Fig. [71 The behaviors can be catego- 
rized into three: (i) the normal bystander effect (NBE), 
i.e. p decreasing with k, (ii) the reverse bystander effect 



(RBE), i.e. p increasing with k, and (iii) the complex 
bystander effect (CBE), i.e. p decreasing and then in- 
creasing with k. With b fixed, for each value of a, there 
appear two critical values of the potential cost c, separat- 
ing the three different BE's. We observe the NBE phase 
for c > CN(a), the CBE phase for CR(a) < c < CN(a), and 
the RBE phase for c < CR(a), respectively. 

To understand how the various BE's appear as the po- 
tential cost c varies, we first consider the simplest case of 
a = in Eq. ^ . For large c, the possibility of successful 
intervention becomes small. This negative effect in the 
link density (small p) is accelerated with k, because the 
/c-dependent interaction term becomes negative. There- 
fore, we expect that p decreases with k (NBE) in the 
stationary state. For small c, the tendency toward suc- 
cessful intervention is accelerated with A: in a positive 
way. So we expect the RBE. 

To be more specific, the dynamics can be simplified by 
the mean field (MF) approximation of Xm at a = 0: 

p{t + At) = 1 + b{k - l)[2p{t) - 1] - c, (25) 

then the stationarity condition yields 

^ , - c 

P-^2 + l~2b{k-l) (2^) 
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with Co = 1/2. One can easily see the NBE phase for c > 
Co, and the RBE phase for c < Co- Numerical simulation 
results are consistent with this MF picture, except that 
Co ~ 0.73. 

The positive a enhances the possibility of selecting the 
witness i connected to the victim (piv = 1) as an inter- 
vener, and hence the fc-hole constraint conies into play. 
The fc-hole constraint, in general, suppresses the link den- 
sity p. So we expect that the NBE phase is not affected, 
because a only accelerates the decay of p with k faster. 

However, the RBE phase may be modified significantly 
in the region where the fc-hole constraint becomes impor- 
tant. Here, the interaction term proportional to b(k — 1) 
competes with a. For sufficiently large fc, we may ignore 
a compared to the interaction term, and then p increases 
with fc. However, for small fc, the effect of positive a 
becomes strong so the fc-hole constraint may make p to 
decrease with fc up to some value of fc. One may expect 
that this non-monotonic behavior of p with fc (CBE) ap- 
pears between the NBE and the RBE phase (Fig. El^d)). 

Note that there is a big jump in the phase bound- 
ary at a = 2b. As discussed in the previous section, 
the fc-hole constraint becomes loose for a < 2b, and be- 
comes strong for a > 2b. As a result, the link density 
decreases discontinuously as one crosses the a = 2b point 
(see, e.g.. Fig. (5]). So it is easily expected that the NBE 
phase expands and the RBE phase shrinks. Moreover, 
the RBE phase effectively disappears for a > 26, because 
it is almost squeezed down to the c = a line where the 
paradise phase (p = 1) starts to appear. 

The complex and reverse bystander effects found in our 
model can support the unusual bystander effects found 
for the nonserious experiments (small c). Moreover, when 
the relation between witnesses is more important than 
the witness- victim relation (a is smaller than &), our 
model predicts that the CBE and the RBE are more com- 
monly observed in nonserious experiments. To confirm 



the robustness of CBE and RBE in reality, more empiri- 
cal/experimental studies are necessary. 



V. SUMMARY 

We have studied how the collective helping behavior 
and bystander effects at the level of society emerge from 
the repeated helper-recipient interactions by means of the 
coevolving helping network model. Its dynamic rules are 
based on the results of social psychological experiments. 
By analyzing the dynamics of the network link density 
(aggregate expectation of help) and its steady states and 
fiuctuations, we present the full phase diagram with var- 
ious active and inactive phases and explore the nature of 
phase transitions between them. Close to the transition 
between the active and paradise phases, the critical be- 
havior turns out to depend on the number of witnesses 
per accident fc. This is due to the link creation suppressed 
by the fc-hole constraint, which also governs the algebraic 
decaying behavior in the paradise phase. 

We have numerically found various kinds (normal, re- 
verse, and complex) of collective bystander effects for 
large fc and proposed the underlying mechanism. The 
normal and complex (non-monotonic) bystander effects 
are consistent with the field study results for the cases 
with serious and nonserious helping situations, respec- 
tively. In addition, we expect the reverse bystander effect 
to occur when the witness-victim relation is sufficiently 
less important than the relation among witnesses for non- 
serious helping situations. 
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